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Abstract: A two-dimensional cellular automaton is introduced to model the 
flow and jamming of vehicular traffic in cities. Each site of the automaton 
represents a crossing where a finite number of cars can wait approaching the 
crossing from each of the four directions. The flow of cars obeys realistic traffic 
rules. We investigate the dependence of the average velocity of cars on the 
global traffic density. At a critical threshold for the density the average velocity 
reduces drastically caused by jamming. For the low density regime we provide 
analytical results which agree with the numerical results. 
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1 Introduction 



Vehicular traffic on highways as well as in cities tends to suffer from a jamming 
transition when the global traffic density exceeds a critical threshold value. The 
phenomena related to traffic jams have attracted the attention of engineers 
and physicists since many years and there exists a large variety of experi- 
mental observations, see fl|dll and references therein. The formation of traffic 
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jams has been investigated by many authors using various methods. As early 
as 1955 Lighthill and Whitham |||| described the spontaneous formation of 
regions of increased car concentrations and shock waves using their approach 
of kinematic waves, i.e. wave motion where the spatial value of the flow is a 
function of the spatial concentration distribution. Prigogine and Herman ]^.|^| 
investigated traffic jams on highways by performing a stability analysis of 
hydrodynamic equations. One-dimensional cellular automata models |J for 
the simulation of traffic flow have been proposed e.g. in [jlO| pT2| . This type of 
models has been studied intensively using massively parallel computers JT3| - 
iTf . Schadschneider and Schreckenberg [[HJ solved a one-dimensional cellular 



automaton model analytically in the mean field approximation, Csahok and 



Vicsek |19| investigated their automaton in the presence of quenched noise. 



Nonlinear wave descriptions have been proposed in |20H22| . Ben-Nairn et. al 
PB| apply a ballistic aggregation process to model the kinetics of clustering 
in one-dimensional traffic flow. Hydrodynamic approaches have been inves- 
tigated by various authors, e.g. [^4],^]. These models are based on certain 
assumptions concerning the "hydrodynamic" properties of traffic flows such 



as the velocity versus density relation P^f27|] , and the viscous terms |2"2] 



One-dimensional traffic has been studied extensively and the developed mod- 
els are sufficiently sophisticated to reproduce experimental observations, i.e. 
the fundamental diagram (throughput versus flow density) of one lane roads. 
Some of the authors (e.g. p8[ ) claim that the size distribution of the jams 
obeys a power law, but very large scale computations have shown that the 
lifetimes of the jams and hence their sizes reveal a characteristic cut-off which 



seems to be no finite size effect of the simulation | l3fl . There are several models 
describing traffic flow in more than one dimension. The analysis of traffic on 
a network with multiple sources and sinks can be found e.g. in [[2~9|-|32"|j . Na- 



gatani |33| describes a cellular automaton representing a two-lane roadway. 



An interesting new approach was currently proposed by Bando et al. |34| who 
express the rules which govern the behavior of the cars (i.e. acceleration and 
deceleration) by a dynamic equation for each car x n = a{V (x n+ \ — x n ) — x n }. 
The interaction of neighboring cars is expressed in terms of the function V 
of the distance between the cars. Using this approach the authors connect in 
some respect the model of the cellular automata with molecular dynamics. 

Two-dimensional cellular automata designed to simulate traffic in a city have 
been investigated by some authors. Biham et al. |!5| proposed a three state 



cellular automaton model where at a given time step each site can be occupied 
by a car moving from South to North, by a car moving from West to East, or 
the site may be empty. The time behavior is ruled by synchronous traffic lights 
at each site, allowing alternatively for vertical or horizontal traffic. Obviously 
there are situations when the traffic jams, i.e. when the traffic light allows 
for driving but the next site is occupied. The model given in p3 and par- 
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ticularly the spatial correlations in the jamming phase have been numerically 
investigated by Tadaki and Kikuchi p6 . A similar model but with "faulty" 
traffic lights was studied in |37[]. Nagatani [3S| investigated the spreading of 
a jam which is induced by an accident using an extremely simple automaton 
rule. Numerically he finds scaling laws for the size of the spreading jam as a 
function of time elapsed since the occurrence of the accident which cannot be 
derived analytically so far. Fukui et al. numerically investigated the evo- 
lution of ensemble averages of the jamming process for the simple automaton 
model described in [fj5[| . For the same model Chau et al. [20] recently gave 
an analytic upper bound (depending on the dimension) for the critical car 
density r] cr , i.e. when the system transits from the "moving phase" into the 
"jamming phase". In two dimensions they found t] cr < 11/12. Obviously this 
upper bound is not close to the value which was observed in numerical simula- 
tions. In the case of Nagatani and Seno there is a set of parallel one-way 
streets, all oriented in the same direction (x-direction) , which is intersected 
by a single perpendicular one-way street in y-direction. 



For fixed car density on the perpendicular street p y they find that the flux 
of cars J x rises linearly with the car density p x until a characteristic critical 
threshold p c x is reached. When further increasing the density p x the flux J x 
approximately remains constant while the average velocity of the cars (v x ) 
drops. This rather sharp change in the overall behavior is due to the forma- 
tion of traffic jams where the flow is irregular and discontinuous. When further 
increasing the density p x the discontinuous character of the flow disappears 
at a second threshold p^ and the flux J x declines linearly. This behavior cor- 
responds to the symmetry of the car density in y-direction and the spaces 
between the cars 1 — p y . Indeed, one finds that the critical densities p c x and 
p x are almost exactly symmetric, e.g. for p y = 0.3 one finds p x . = 0.31 and 
I'', = 0-69. 



Cuesta et al. [0,0 have been the first who reported on simulations of a two- 
dimensional automaton where the cars are allowed to change their direction. In 
their model each crossing can be occupied by one car or it can be empty. Each 
car is assigned a preferred direction; the parameter Wi(r) is the probability 
that the i-th car at the position r moves in the next time step horizontally 
and 1— Wi(r) is the probability to move in vertical direction. Horizontal motion 
is allowed at even time steps, vertical motion at odd times. There have been 
defined two variants, A: there are only one-way streets directed from South to 
North and from East to West. Half of the cars is given the trend Wi(r) = 7 and 
the other half is given Wi(r) = 1 — 7. In the second variant B there are one- 
way streets pointing alternatively South — > North and North — > South, and 
West — > East and East — > West, respectively. The trends of the cars which 
are subdivided into four equal groups, point into one of the four possible 
directions. 



3 



We should mention that several authors claim that at least the one-dimensional 
traffic flow problem is closely related to sand flowing in pipes, e.g. 48 . 



In the present paper we are concerned with a two-dimensional cellular au- 
tomaton designed to simulate traffic in a city. As a special feature of our 
model we emphasize the fact that we use rather realistic traffic rules which 
will cause a slowing down of the average traffic velocity and finally will lead 
to the collapse. 



2 The model 



In our model the city is represented by a set of L streets in horizontal direction 
crossing L streets in vertical direction. Cars are allowed to move in both 
directions, i.e. we have no restrictions to one-way streets. The crossings define 
a two-dimensional cellular automaton, each of the L 2 crossings is represented 
by a site of the automaton. We assume periodic boundary conditions in both 
directions. Fig. [I] shows a schematic plot of an automaton site (a crossing). At 
each crossing there are four queues of maximum length Q filled by cars coming 
from one of the four directions, representing the finite space where cars can 
move freely between crossings in realistic urban traffic systems. Hence, at each 
crossing there are allowed at most 4 Q cars. Each of the cars has a desired 
direction, i.e. to the right, to the left or straight on, which has to be chosen 
according to a certain rule. We will discuss this point in detail below. Within 
each time step the first of the cars in each queue can go to the next crossing 
provided that it does not have to give way to one of the other cars at the same 
site. In case it has to give way to another car it will not move during this time 
step. In our model we assume the simple and realistic rule that a car has to 
give way in the case that another car occupies the same crossing at its right 
hand side. Moreover a car has to stop if it intends to turn left and there is a 
car at the opposite side of the crossing which goes straight on or turns to the 
right. If all of the four top positions of a crossing are occupied by four cars 
one of the cars will be selected randomly and will then be allowed to move in 
the current time step. Insofar we have chosen realistic rules which are current 
law in many countries. 

Each of the N cars in the system starts at a randomly selected crossing, 
its desired direction (left, straight on, right) will be determined according to 
some rule (discussed below). Within each time step a car maximally moves a 
distance of one lattice cell, i.e. it tries to transit to one of the next neighbors 
of the current site provided the following three conditions hold: 

(i) It does not have to give way. 

(ii) It is placed at the top position of its queue. 
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Fig. 1. The schematic plot of a crossing 
^iii) There are less than Q cars standing in the destination queue. 



Otherwise it rests. Obviously, if there is only one car in town it never has to 
give way nor will it find another car in front of itself. Therefore the average 
velocity^ (v) = 1 in units of lattice cells per time steps. 

As mentioned above we have to discuss the rule which determines the desired 
direction within each time step. There are at least three reasonable rules: 

(i) The desired direction for each car is selected at random within each 
time step not regarding whether the car moved in the previous time step 
or stopped. This rule is rather artificial since provided a jam occurs it 
can dissolve due to the fact that the cars might give up their previously 
desired direction when they cannot move. 

(ii) Each car is assigned a destination site on the lattice. Once a car reached 
its destination it will be assigned a new randomly chosen destination. 
This rule is problematical since one has to define the detailed path a car 
follows to reach its destination. A natural choice would be to chose a path 
so that the car has to change its direction only once on its way between 
its starting position and its destination. In this case, however, the overall 
behavior of the system would depend on the size of the lattice L since one 
can check that the probability to be allowed to move straight on is higher 
than to be allowed to turn. Cars which intend to turn either left or right 
(on the average) have to give way more frequently than those which go 
straight on. The larger the system the smaller becomes the probability 
per time step to turn. For the case L — > oo the probability for a car in a 
given time step to turn approaches zero. Hence this rule does not allow 
for size-independent results, and the results of a simulation using this 
rule apply only for one fixed lattice size L. 

(hi) Each car is assigned a desired direction during the initialization. But 

1 Here and in the following we describe the ensemble average of a value x by the 
symbol (x), whereas the time average by x. 



5 



now, only after moving one step to the aimed destination a new desired 
direction will be chosen at random, i.e. the primarily chosen desired di- 
rection is maintained by each car until all of the three conditions hold 
which allow the car to move. This rule has several advantages compared 
with the other rules. First we avoid the size dependent behavior of the 
previous rule, and second we do not have long range correlations between 
the lattice sites for the case of low car density, which would be caused by 
the second rule. Hence the behavior of the system can be described by a 
Markov process, which will be the basic assumption of our analytical cal- 
culation in section |j. An analytical description which takes into account 
the long range correlation between the sites coming from the rules which 
determine the path of the cars seems to be very difficult, except perhaps 
for the (trivial) case L — > oo, when the probability of the cars to turn, 
and hence the traffic rules connected with the turn, can be neglected. In 
the following we always refer to the third rule. 



3 Numerical results 

In this section we present the results found by numerical simulations of a 
cellular automaton which represents the traffic rules described above. 

The two-dimensional automaton consists of L x L sites on a rectangular lat- 
tice. To eliminate boundary induced effects we have chosen periodic boundary 
conditions in both directions, i.e. we assume the topology of a torus. We sim- 
ulated systems of size L = 20, 30, 50, 100 and no striking difference in the 
behavior of the system was observed. Hence, we exclude finite size effects. The 
fact that already comparatively small lattices exhibit this independence from 
their size seems to indicate the existence of only weak spatial correlations be- 
tween sites separated by a large distance. For high densities, i.e. close to or 
above the critical density r] cr , this statement does not remain valid. In fact, we 
expect the existence of long range correlations. However, in the present paper 
we are mostly concerned with the low density region and we will not devote 
ourselves extensively to the behavior of the automaton for high car density. 

Each of the L x L automaton sites had its own four queues of length Q = 10 
connecting the site with one of its four neighboring sites as described in the 
previous section. For each value of the density i] = where N is the number 
of cars, we started the system at randomly chosen initial conditions. Each car 
was assigned the initial site, the queue in this site (i.e. whether it comes from 
North, East, South or West) and the desired direction (left, straight, right) for 
the next step. When a car transits to an empty site from a certain direction 
it occupies the first empty place in the queue. If it is at the top position of 
the queue it can transit to the next site or it has to wait, depending on the 
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cars coming from the other directions of the same site, and depending on the 
traffic rules which apply to the given situation. If a car transits to another site 
and if there are other cars located in the previously mentioned queue, all of 
them advance by one position in this queue hence, moving up "closer to the 
crossing" . 

After starting the simulation for a given number of cars, i.e. for a given density 
T), we evolved the system 10,000 time steps without recording the statistics in 
order to let the system forget about the initial preparation. Typically a "sta- 
tionary" behavior was found after only a few hundred time steps. In density 
regions very close to the transition value however, the relaxation of the system 
lasts longer which hints at some phase-transition-like phenomenon. The mean 
velocity of the cars was then computed using 2,000 time steps. We varied the 
car density beginning from 77 = to rj = 3.8 in steps At] = 0.01. 

The results of the simulations are shown in fig. |^. Like in simulations of cellular 
automata performed by other authors with different traffic rules we found a 
smoothly decaying function in the low density region. At a certain threshold 
value r] cr fa 1.6 the character of the traffic flow changes abruptly and the 
system transits into the clogged regime. But different from other systems we 
find a remnant velocity in the clogged regime. This is due to the rule that 
in case all four of the queues are occupied there will be one car selected 
randomly which then will be allowed to move. Therefore, in our model there 
exist situations where a jamming at a crossing with four queues occupied may 
still dissolve. 

The curve in the small density region is very smooth, i.e. there are almost 
no fluctuations. Hence, we can assume that long range correlations do not 
play a major role there. In the following section |3] we will present a statistical 
description appropriate to model the flow in the low density regime, i.e. for 
V *C i] cr , which was inspired by this observation. 



4 A statistical description of vehicular traffic 

Let the number of cars be denoted by N and the number of cells respectively 
crossings building the torus by M = L x L where L is the common perimeter 
of the torus. For the global density we choose the symbol 77 = N/M. 

The quantity we focus our attention on is a time averaged^ and normalized 

2 When performing the simulation the parameter T is chosen sufficiently large in 
order to scan the state space respectively the limit set with sufficient accuracy. 
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Fig. 2. The simulation result: The mean velocity v s i m of N cars on a periodic lattice 
of size L x L vs the global traffic density rj = N/ (L x L), Except for values around 
the critical density rj cr the curve is rather smooth. In the vicinity of the transition 
zone critical fluctuations, i.e. temporary jams, cause an irregular relationship. This 
can be understood from the fact that the typical lifetime of jams becomes of the 
same order as the simulation time. 
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where A(t) means the total number of cars prevented from moving by traffic 
rules at time t which depends on the exact situation experienced during the 
simulation. 



For a description of the dynamics in the low density region f| C 1 we apply 
the following arguments: 

(i) An effective description is possible which separates the process of clus- 
tering at the M possible crossings from the process of obeying the traffic 
rules. Due to this independence assumption the time average ([!]) ex- 
tracted from a simulation can be decomposed into two averages: first a 
time and cluster average related to the traffic rules and second a time 
average with respect to the clustering process. 

(ii) Given i cars meeting at one crossing a variety of situations is possible. 
Every car is approaching the crossing from either North, East, South or 
West and intends to go straight, to turn left or to turn right. In the long 
run all possible situations will have equal probability. Hence, the average 
number of cars which must wait, denoted by 5i with < 5{ < i, can be 
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computed by combinatorial reasoning. The result of this computation is 



collected in table |A.1| for i = 1, . . . , 8. A sketch of the derivation is given 
in appendix [A]. 

For clusters of size i > 8 the 8i are not shown since their contributions 
are suppressed by very small cluster probabilities, 
(iii) Central to our description is the identification of cluster distributions, 
i.e. we write k = (ho, ki, . . . , /cjv) meaning that ki denotes the number 
of cells in the system where i cars are meeting. Of course, there are two 
boundary conditions, namely 

N 

E^ = M (2) 

and 

N N 

'£ik i = '£ik i = N . (3) 

i=0 i=l 

We now assume that for the low density regime rj <C 1 the process of 
clustering is ergodic which means for almost all initial configurations the 
time average can equivalently be represented by an average according to 
an invariant probability distribution. Note that this probability distri- 
bution is defined over the space of all configurations k compatible with 
the boundary conditions (|2|) and ([$[). In order to derive this distribution 
we additionally assume that the dynamics effectively is equivalent to a 
process of independently distributing A" cars (balls) among M cells with 
crossings (urns) according to the equidistribution. This means we assume 
the discrete dynamics to be a Bernoulli process. Then the probability dis- 
tribution reads 

M' AT' 

Pit) = jf • (4) 

M N II kj\] 

3=0 

To understand formula ([§) in detail one has to realize that there exist 
Ml I [fco'^i! • • • ktf\] different ways to index the cells with numbers 0, 1, . . . , A" 
since all cells containing the same number % of cars are not distinguish- 
able. Then there exist N\/ (0!) fe °(l!) fcl . . . {N\f N ways to index the cars 
with numbers 1, . . . ,N taking care of the fact that cars within an identi- 
cal cell are not distinguishable. Finally the equidistribution with respect 
to all configurations yields the factor M~ N . This result was obtained by 
von Mises (1939) @. 

Now the time average (A)j is replaced by the ensemble average 

* N 
k i=l 

The sum with respect to i reflects the average delay due to the traffic rules 
at the crossings and the sum with respect to all configurations obeying (0) 
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and (H) - which is indicated by the symbol * - accounts for the cluster statis- 
tics. The factor a, accounts for the fact that the dynamics of the real process 
differs from a Bernoulli-process. In general there may exist spatial as well as 
temporal correlations. The spatial correlations become important for densely 
filled lattices. The dynamics which control the traffic within one cell is strongly 
influenced by the cars moving in neighboring cells. Finally the emergence of a 
traffic jam clearly expresses these correlations. On the contrary temporal cor- 
relations play the dominant role for very sparse system, i.e. in the low density 
region. Consecutive situations are not completely statistically independent. 
Let us briefly explain this effect considering only cluster of size i = 2: When 
the cars are randomly distributed amongst the lattice sites there is a certain 
probability that two cars p and q occupy the same site and hence form a clus- 
ter of size i = 2. Assuming statistical independence this probability is time 
invariant but for the cellular automaton the situation is quite different. 

If two cars meet at a site there exist two possibilities: either both cars continue, 
or one of them continues and the other one has to give way. There is no chance 
that both cars meet again one time step later. In the extreme case of only two 
cars, i.e. N = 2, the probability for them to meet by random distribution is 
thus effectively reduced by a factor of two in comparison with the automaton 
dynamics. Given only a small number A" of cars, i.e rj = N/M <C 1, and 
provided only two of them p and q meet occasionally one has to consider in 
the following time step a system where only N — 2 cars feel the interaction 
with (in our case: can meet occasionally) N — 1 other cars and 2 cars which 
can meet only N — 2 cars, i.e. p is not allowed to meet q and q cannot meet 
p. Hence, the system has some kind of memory and therefore it is not ideally 
Bernoullian. 

Obviously this effect plays a crucial role only for small A" since otherwise the 
relative difference between A^ — 1 or A^ — 2 is negligible. Speaking in terms of 
density this means that the effect is substantial only for very small densities. 
Then however, clusters larger than i = 3 are suppressed by small likelihood. 
The net effect of these correlations can be taken into account through counting 
all possible configurations which can be realized after two cars have met. This 
procedure is similar to the method employed for the determination of the Si (as 
explained in appendix [A]). Neglecting contributions related to clusters of size 
larger than three, which is safe for very small densities (see fig |3|, explanation 
below) means ai = 1 for i > 4. From simple geometric considerations we find 
cui = 1, c*2 = 1/2 and «3 = 7/4. 

Inserting expression eq.(|) into ([!]) yields 
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Fig. 3. The average number of clusters of size i normalized to the lattice size 
M = L x L as a function of the global traffic density r/ = N/M. The lines represent 
eq. (9) while the points correspond to numerical data achieved by a simulation on a 
50 x 50 lattice (torus). The curve i = gives the probability for a site to be empty. 

and performing the sum over all kj (j ^ i) we find 



i=l JV ki=0 

N ~r 

i-E^(if ( ). 



i=l 



AT 



(7) 



The distribution p(ki) - declared for i = 0, 1, ... , [N/i\ - can be calculated 



using the inclusion-exclusion principle [p0| , |51| . This derivation is performed in 
appendix [B] and the result reads 

i\n m / ,\ (m — jY N -rt 



j=k 



We see from (0) that the ingredients we actually need are the first moments 
(Ki) of the cluster distribution. They can be derived calculating a generating 
function for the (descending factorial) moments, denoted h\ m \z,x) pOfl . This 
calculation is performed in appendix [Cj The result reads 



/jV\ 1 / 1 

<*<> = (;) P^('-m) < 9 > 
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In fig. m we have depicted the expectation value of the cluster distribution (Ki) 
(eq. divided by the system size M = L x L. The solid lines correspond to 
the function given in eq. (|J) and the points represent related data taken from 
our simulation on a 50 x 50 torus. The values for % = give the probability for a 
site to be empty. For i = and i — 1 the analytical and numerical results agree 
perfectly, for larger i the discrepancies between the solid lines and points are a 
direct consequence of the deviations from the assumed independent behavior 
as explained above. Note that the probabilities in fig. ^| are plotted using a 
logarithmic scale. The observed discrepancies for clusters of size larger than 
3 will hardly have any influence on the overall behavior of the automaton. 
For values rj > r] cr the simulation data illustrate a breakdown for the small 
occupation numbers (except % = 0) due to the emergence of jams which are 
clusters of high order. At the transition point the occupation rates abruptly 
drop by three orders of magnitude. The precise position of this transition 
point is masked by the finite size of N used in the simulation. Increasing 
the car density rj beyond the blurred critical zone results in a slow rise of the 
occupation number for the small clusters and simultaneously, in a slow decline 
of the number of empty sites. This is due to the increasing number of cars still 
moving around the jammed crossings. 

Inserting equation (||) into (|7|) yields the following formula 

v = l-Y ^2 ( iV ) _i_ (1- — ) . (10) 

The last step will be to substitute in (|i0|) the density r] for the number of cars 
according to N = rjM which results in 
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The main result of our statistical description, namely the mean velocity as 
a function of the car density (eq. (|TTD), is plotted in fig. [|. The solid line 
shows the function given by eq. (|TT|) and the points are the data taken from 
the simulation on the 50 x 50 torus (compare fig|| In the low density region 
we find a nice agreement. The closer one approaches the transition point the 
more eq. (JTTJ) overestimates the simulation results. This can be understood 
by recalling the fact that for densities close to the critical value the effective 
decomposition of the automaton dynamics into the cluster dynamics and the 
single site dynamics becomes inappropriate due to formation and dissolution 
of short lived jams (like critical fluctuations). They tend to increase the oc- 
cupation number of larger clusters and hence, lead to an effective decrease of 
the average velocity. 
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Fig. 4. The average velocity vs the global traffic density for the low density regime 
(non jamming traffic). The curve displays the result from equation (11), the points 
show data taken from the simulation. 

In the thermodynamic limit - which means M, N — > oo and keeping n = N/M 
constant - we drop terms of order 0(M~ l ) and end up with 

r]M J- 

v « 1 — — — i] ( > exp(— 77) (12) 

i=l % - 




where the exponential corresponds to the dominant contribution of the last 
factor in (ITT]). This formula can be expressed by a series expansion in powers 
of 7] 

l) k 8~a, (f _ 1+Jfe) 

_ 5 3 a 3 \ 2 _ / 5 2 a 2 _ 5 3 a 3 5 4 a 4 
6 r V 4 6 24 





Note that there is no constant term in the double sum since 5\ = 0. Equa- 
tion ( |1~3| ) can be truncated at a given order of i] hence, giving rise to an 
approximation scheme. In fig. || we plotted the linear (dashed), the quadratic 
(dotted) and cubic (long dashed) approximations together with the full for- 
mula (fat solid). The range of reliability visibly increases with increasing the 
order of approximation. On the other hand this scheme clearly illustrates that 
the functional relation between mean velocity and global density definitely is 
nonlinear, except for very small car densities. 
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Fig. 5. The average velocity as derived from the statistical description (eq. (9)) (fat 
line) and truncated series expansions: linear (dashed), quadratic (dotted) and cubic 
(long dashed). This plot clearly illustrates a nonlinear relationship between average 
velocity and global density. 

5 Conclusion 

We investigated the behavior of a two-dimensional cellular automaton with 
periodic boundary conditions to simulate traffic flow in cities. The automa- 
ton mimics realistic traffic rules which apply to our everyday experience of 
vehicular traffic. In agreement with similar models we found a slow decay for 
the mean velocity (v) as a function of the global traffic density rj. At a crit- 
ical threshold value rj cr the mean velocity collapses abruptly and the system 
transits into another regime of global behavior which we call the jamming 
regime. Beyond the critical density rj cr the average velocity is very small and 
it declines further when increasing the density. By simulating different sizes 
of automata we could exclude the influence of finite size effect provided the 
density value is not rather close to the critical density. 

Applying combinatorics and statistical methods for the description of the sys- 
tem in the low density regime the analytical calculation performed in section |] 
yielded the average car velocity as a function of the car density - (v) (77) - which 
nicely agrees with the values of the numerical simulations in section [|. Since 
the derived functional relationship between mean velocity and global traffic 
density was based on very general assumptions (conservation of the number 
of cars, weak spatial correlations in the low density regime, ergodicity) we ex- 
pect this description to be valid for a wider class of traffic systems. Note that 
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the specific structure of traffic rules enter the description only through the 
average delays St. The only nonlocal ingredient arises from the restriction that 
cars have to stop in case the next desired queue at a neighboring site is totally 
filled up. But this nonlocal character only becomes substantial for densities 
close to the critical value where it causes long range spatial correlations. 



Several authors (e.g. [p2jD assert that the average velocity v in the low den- 
sity regime is a linear function of the density rj and indeed the simulation 
results seem to support this observation. A more detailed analysis however 
reveals that this function definitely is nonlinear. An analytic expansion shows 
(eq. fll3|) , fig. HD that the results become dramatically wrong when truncating 
the formula after the linear term. Surprisingly the analytic description gives 
good results even for densities not too small where the nonlocal effects of the 
dynamics and hence, long range correlations spoil the basic assumptions of 
our treatment. 

For still higher values of the density 77 the results of our simulation agree with 
results reported in literature (e.g. |52| and many others), but no longer with 



our approximate theory. Above a critical density we observe an abrupt transi- 
tion of the system into the jammed regime where the averaged velocity is close 
to zero. This regime however is beyond the scope of our statistical description 
and has to be investigated starting from other approaches (e.g. nucleation 
processes). 
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A Calculation of the average number of waiting cars Si for clusters 
of size i = 1 and % — 2 



In table [A.l| we present the average number Si of cars which are prevented 



from moving given i cars meeting at a crossing. Obviously the number Si is 
confined to the interval [0, i] . 

The average is performed with respect to all possible configurations which 
are assumed to possess equal probability. The case i = 1 is rather simple: If 
there is only one car at the crossing it has never to give way hence, Si = 0. 
The case i = 2 is not that trivial and the computation of S2 requires some 
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Table A.l 

The average number of cars which are stopped <5j when i cars meet at a crossing. 

#cars meeting i 5i 



1 

1 


n n 


2 


0.638889 


3 


1.479167 


4 


2.467014 


5 


3.524740 


6 


4.604709 


7 


5.683838 


8 


6.752964 



combinatorial effort. Since each of the cars can wait at one of the four sides of 
the crossing and can proceed in one of the three directions (left, straight on, 
right) we have 4x3x4x3 = 144 different situations of equal probability. 
Each of these situation can be assigned to one of the situations in the left 
column of table |A.2|. We now explain the columns of table [A.2|: 



(i) A sketch of the configuration given by a related graph. 

(ii) The number of the different realizations that relate to the graph in the 
first column assuming first that the cars are distinguishable. The first 
factor A in the form A x B x C originates from the symmetry according 
to rotation by ~7r, n and ~7r. The second factor B denotes the number of 
choices for the first car (which is now assumed to come from South) and 
the last factor C gives the number of situations possible for the second 
car. 

(iii) When we take into account that the cars in fact cannot be distinguished 
the number of different events that belong to the figure in the first column 
has to be multiplied by the factor given in the third column. 

(iv) The total number of events for the whole class denoted by the graph, 
i.e. (column 2x column 3). 

(v) The number of cars stopped by the traffic rules at this crossing in one 
of the possible realizations (either one or none). 

(vi) The total number of cars which have to stop related to the whole class 
denoted by the graph, i.e. (column 4x column 5). 



From the last line in table [A.2| one sees that there are 144 different situations 
of equal probability which amounts to 288 cars. Furthermore we see that 92 
of those 288 cars have to stop. Hence, the average number of the cars stopped 
for 2-clusters is T 2 = 92/144 « 0.639. 
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Table A.2 

All situations which might occur when two cars meet at a crossing, their frequency, 
their symmetry according to rotation and due to indistinguishability of the cars and 
the number of cars which have to stop in the current situation (explanation in the 
text). 

graph rot. symm. dist. ^events waiting ^stopping 



4x3x3 1 36 



4x1x2 2 16 



4x1x1 2 



4 x 1 x 3 2 24 



4 x 1 x 3 2 24 



2x1x1 2 



4x1x2 2 16 



2x1x1 2 



4x1x1 



2x1x1 



36 



16 



24 



16 



total 



144 



92 



For the clusters of size i = 3 the average number of cars which are prevented 
from moving <5 3 can be calculated in analogy however the corresponding table 
contains about 8 times as much situations as table |A.2| . Therefore we do not 
want to present it here. The result is 5% = 1.479. For cluster of size larger than 
3 the corresponding 5i have been calculated by the method of complete enu- 
meration using a computer. These empirical results are collected in table |A.l 
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For % = 1 . . . 3 they are identical with the analytical results. 

For completeness sake we remark that for large cluster sizes [i > 8) there 
exists a negligible probability that one of the directions of the crossing is not 
occupied. Typically there are many cars waiting in each of the four queues. In 
those very likely cases our traffic rules allow exactly one (randomly chosen) 

cars to drive, all other i — 1 cars have to stop. Therefore we find 5i i — 1. 
Certainly the contributions of such big clusters will only play a minor role in 
our calculations due to very small likelihood. 



B The derivation of the cluster probabilities 



The probability p (ki) to find exactly k sites which are occupied each by i 
(independently moving) cars in a system of M crossings and iV cars is given 
in eq. (§). The derivation of this formula employing the inclusion-exclusion 
principle will be performed in this appendix. Note that the problem to find 
p (ki) is different from the trivial problem to find the probability for the ki 
crossings which are occupied by at least i cars. 

The inclusion-exclusion principle relates the probabilities for a finite number 
of sections of events to the probabilities of an exact number of events. 



Let (Q, P) be a probability space, Ai,...,A^ be events, and for arbitrary 
{mi, . . . ,rrij} C {1, . . . , iV} let P(A mi n ... Pi A m . ) be known probabilities. We 
define 

Sj= E P(A mi n...nA mj ) (B.i) 

mi,...,m,j 

Moreover, let be B n = {uj 6 Q : uj G A m for exactly n values of k}. Then the 
inclusion-exclusion principle asserts that 

P(Bn) = E (-l)&'- n > ( J )Sj (B.2) 



n 



Applied to our case the Sj are the probabilities that j cells each contain i cars 
and the rest of the cars, i.e. (N — ji) cars, are distributed arbitrarily among 
the remaining (M — j) cells; this probability can be derived with ease and 

reaC ^ S " ■ ( N ") 

s > = (i\y(N-ji)\ (m) i 1 ~ m) (R3) 

To explain this formula we first index each of the N cars with the numbers 
1, . . . ,N which yields the factor N\. Since the % cars within each of the j cells 
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are not distinguishable we have to divide this factor by (i\y and since the 
remaining (N — ji) cars are not distinguishable too additionally by (N — 
After having numbered the cars we subsequently fill the first cell with cars 
1, . . . , i, the second cell with cars (i + 1), . . . , (2i) and so on. In this way we 
distribute the cars numbers 1, . . . , ji among the urns 1, . . . , j yielding the fac- 
tor M(~ Jl \ The remaining (N —ji) cars are distributed successively among the 
remaining M —j cells at random which explains the factor [(M — j) j 'My~^\ 



Insertion of this probability into the inclusion-exclusion formula ( |B.2|) yields 

M 

E 

j=ki 



P(h) = E (-1)^) (£) Si 

ht K ] UJ (N-ji)\ 1 ] 



Because of the generalized definition of factorials - using the gamma function 
- the sum effectively only ranges from ki up to \_Mji\ which is sensible. 



C Calculation of the first moments of the cluster distribution (Ki) 



In the following we will derive the mean value of the cluster distribution (Ki) 
which is given in eq. (||). 

The generating function for the (descending) factorial moments of the cluster 
distribution we use is 



H. 



{M) t 



z,x) 



E E ^x ki p(h,N) 

N=0 ki=0 



(C.l) 



where p(ki, N) is the probability to find a value ki for the stochastic variable Ki 
when trying with iV balls. The benefit of such a complicated looking generating 
function is that the sums can be performed yielding an analytical expression 
namely 



n M 



exp H — -{x 



(C.2) 



For the explicit derivation the reader is referred to the book by Johnson and 
Kotz 3 p. 116ff |50| . Moreover, this function is related to the (descending) 
factorial moments of the stochastic variable Ki according to 



' d r 
dx 7 



H 



x, z) 



x=l 



one has to identify j = i,m = M, n, = N, Mi = Ki,g = ki, Pr[Mj = g\n] = p(ki) 
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"£ki(ki-l) ...(h-r+^pik^N) 

N=0 k t =0 



= E*" - 1) • • • - r + 1) 

AT=0 



(C.3) 



hence, the first moment (.£Q) is given as the coefficient of z N in 



AH d 



M N dx 



x=l 



(C.4) 



Consequently we calculate 
N\ ( d 



M N \dx 



exp 2 +— (x - 1) 



iV! 



x=l 



exp 1 



j(Af-l)] 



?,! 



- AH_ (M-l) fc (fc+i 
^ i! fc! M^- 1 ) 



(C.5) 



and finally arrive at 

N\ (M-l)^"*) /AT 



i!(JV-i)! M^- 1 ) 



1 v(JV-<) 



(C.6) 
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